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Abstract. Typical heat engines exhibit a kind of homotypy: The heat exchanges 
between a cyclic heat engine and its two heat reservoirs abide by the same function 
type; the forward and backward flows of an autonomous heat engine also conform to the 
same function type. This homotypy mathematically reflects in the existence of hidden 
symmetries for heat engines. The heat exchanges between the cyclic heat engine and 
its two reservoirs are dual under the joint transformation of parity inversion and time- 
reversal operation. Similarly, the forward and backward flows in the autonomous heat 
engine are also dual under the parity inversion. With the consideration of these hidden 
symmetries, we derive a generic nonlinear constitutive relation up to the quadratic 
order for tight-coupling cyclic heat engines and that for tight-coupling autonomous 
heat engines, respectively. 
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1. Introduction 

As an important sonrce of power, heat engines are crncial to onr hnman activities. 
It is necessary to investigate their energetics in onr times of resonrce shortage. The 
classical eqnilibrinm thermodynamics provides a powerfnl tool to investigate the ideal 
heat engines consisting of reversible processes. The concept of Cannot efficiency is one 
of the cornerstones of thermodynamics. It serves as the npper bonnd for efficiencies of 
heat engines working between two heat reservoirs. 

Cannot efficiency is achieved only for heat engines operating in eqnilibrinm or 
qnasi-static states, which operate infinitely slow and yield vanishing power ontpnt. 
The practical heat engines usnally complete thermodynamic cycles in finite periods 
or operating at finite net rates, which give rise to the development of finite-time 
thermodynamics. The most elegant resnlt on this topic is the efficiency at maximnm 
power for the Cnrzon-Alhborn endoreversible heat engine m. which reads rjcA = 
1 — ^Tc/Th with Th and being the temperatnres of the hot reservoir and the 
cold one, respectively. This resnlt has attracted mnch attention from physicists and 
engineers [4ll37]. The previons researches reveal that the Cnrzon-Alhborn efficiency 
(Vca) is also recovered, or at least approximately recovered, in a lot of heat engines snch 
as the stochastic heat engine |30], the Feynman ratchet [31], the single-level qnantnm 
dot engine [32], and the symmetric low-dissipation heat engine [33] . 

The connection between finite-time thermodynamics and linear irreversible 
thermodynamics was proposed by Van den Broeck [M]. The constitntive relation, 
which is defined as the relation between the generalized thermodynamic ffnxes and 
forces, is one of the central formnlas in irreversible thermodynamics. The constitntive 
relation is assnmed to be linear and restricted by the Onsager reciprocal relation [3H] 
in linear irreversible thermodynamics. Using the linear constitntive relation. Van den 
Broeck fonnd that the efficiency at maximnm power for tight-conpling heat engines 
is half of the Carnot efficiency. However, we observed that many heat engines, snch 
as the Cnrzon-Alhborn endoreversible heat engine [3], the Feynman ratchet as a heat 
engine [31] and the single-level qnantnm dot heat engine [32], can not be precisely 
described by the constitntive relation for linear response. There exist higher order 
terms in the relations between generalized thermodynamic ffnxes and forces for these 
heat engines. This conflict might lead to inappropriate resnlts when we try to investigate 
the energetics of heat engines in a higher precision snch as the nniversal efficiency at 
maximnm power np to the qnadratic order [351137] . Thns, it is an nrgent task for ns to 
seek a generic expression of nonlinear constitntive relations for heat engines. 

On the other hand, based on onr early researches we noticed that typical heat 
engines exhibit a kind of homotypy, which had not been tonched in previons literatnre. 
We observed that the heat exchanges between a cyclic heat engine and its two 
reservoirs abide by the same fnnction type; and that the forward and backward 
flows for an antonomons heat engine also conform to the same fnnction type. This 
homotypy mathematically reflects in the existence of hidden symmetries: The heat 
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exchanges between the cyclic heat engine and its two reservoirs are dnal nnder the joint 
transformation of parity inversion and time-reversal operation [see PT-symmetry flT^ 
in Sec. 13.Ij : the forward and backward flows in the antonomons heat engine are also 
dnal nnder the parity inversion [see P-symmetry fl3^ in Sec. 14.1] . It is still nnclear 
what constraints will these hidden symmetries impose on the expressions of nonlinear 
constitntive relations for heat engines. 

In this work, we hrst revisit the generic model for tight-conpling heat engines 
proposed in onr previons work [Ml- Then we will investigate the hidden symmetries 
in heat engines, based on which we derive the nonlinear constitntive relations np to 
the qnadratic order for cyclic heat engines and antonomons heat engines, respectively. 
These results are also conhrmed by typical models of heat engines. 

2. Generic model 

In a conventional setup of heat engine, the engine absorbs heat Qh from a hot reservoir 
at temperature Th and releases heat Qc into a cold reservoir at temperature Tc per 
unit time. Simultaneously, it outputs power W against an external load, which may be 
further expressed as 

W = Qh-Qc, (1) 

according to the energy conservation. 

For the engine operating in a hnite period or at hnite rate rather than in a quasi¬ 
static state, the contribution of the interactions between the engine and the reservoirs 
should not be ignored j36]. We introduce two nonnegative weighted parameters Sh and 
Sc satisfying -t-Sc = 1 to represent the degree of asymmetry of the interactions between 
the engine and the reservoirs. With these parameters, the weighted thermal flux Jt and 
the weighted reciprocal of temperature [3 may be dehned as 

Jt — ^hQc T ^cQhy (2) 

and 

(3 = Sh/Th 3- Sc/Tc, (3) 

respectively. The values of Sh and Sc depend on specihc models. In particular, 
Sh = Sc = 1/2 indicates a special situation in which the engine is symmetrically interacts 
with the two reservoirs. 

From ([1]) and ([2]), we obtain the heat fluxes 

Qh = Jt + ShW, and Qc = Jt - ScW. (4) 

Based on these relations, we obtain a rehned generic model as shown in Fig. [H In this 
new physical picture, the engine absorbs heat Qh from the hot reservoir per unit time 
while an amount of heat ShW is converted into work due to the interaction between 
the engine and the hot reservoir. A thermal flux A flows through the engine. Then 
the engine releases heat Qc into the cold reservoir per unit time while an amount of 
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heat ScW is converted into work due to the interaction between the engine and the cold 
reservoir. The influence of the relative strength of the interactions between the engine 
and the reservoirs is explicitly included in this picture. 



Figure 1: Refined generic model of a tight-coupling heat engine (reproduced according 

to [36]). 

The generalized thermal force conjugated to Jt may be defined as 

X* = 1/T, - I/Th. (5) 

The definitions of generalized mechanical flux Jm and mechanical force Xm depend on 
the type of heat engines. For a cyclic heat engine, they may be defined as 

Jm = 1/to and Xm = -/dW, (6) 

respectively, where to is the period for completing a whole thermodynamic cycle. It 
should be noticed that the sign of to is of physical meaning, to takes a positive sign 
when the thermodynamic cycle operates as a heat engine, while it takes a negative sign 
when the thermodynamics cycle operates as an refrigerator. W is the work output in 
each cycle. For an autonomous heat engine operating in a steady state, Jm and Xm may 
be defined as 

Jm = r and Xm = -/dw, (7) 

respectively, where r is the net flow and w denotes the elementary work in each 
mechanical step. 

With definitions (EJ-jITj), it is easy to verify that the entropy production rate of the 
whole system may be written in a canonical form 

O' = Qc/Tc ~ Qh/Th = JmXm + JtXt- (8) 

In this paper, we focus on a tight-coupling heat engine, in which the heat-leakage 
vanishes so that the thermal flux is proportional to the mechanical flux: 

Jt = iJm (9) 

where ^ is named coupling strength. Its physical meaning may be interpreted as the 
elementary thermal energy flowing through a cyclic engine in each period or that flowing 
through an autonomous heat engine in each spatial step. Define affinity 


A = Xm + iX^, 


( 10 ) 
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which represents the competition between the generalized mechanical force and the 
generalized thermal force. Then, the entropy production rate ([8]) may be further 
expressed as 

a = JmA. ( 11 ) 

3. Cyclic heat engines 

In this section, we will investigate the hidden symmetry in cyclic heat engines and its 
influence on the constitutive relation for nonlinear response. 

3.1. VT-symmetry for cyclic heat engines 

By analyzing typical models of cyclic heat engines in the literature such as the low- 
dissipation heat engine |33], the Curzon-Alhborn heat engine [3], and the revised Curzon- 
Alhborn heat engine , we observed that these models exhibit a kind of homotypy—the 
expression of the heat absorbed from the hot reservoir and that released into the cold 
reservoir abide by the same function type. For example, in low-dissipation heat engine, 
the heat absorbed from the hot reservoir and that released into the cold reservoir can be 
expressed as Qh = Th{AS — S/j/t/j) and —Qc = Tc{—AS — Hc/h), respectively, which 
are exactly of the same function type. Here we have not explained the meanings of 
physical quantities in these two equations. The details can be found in Ref. |33]. The 
other examples such as the Curzon-Alhborn heat engine and the revised one are fully 
discussed in Secs. 13.31 and l3^ We remind the reader to note equations ([23]), fl2Tl) . and 

m- 

This homotypy mathematically reflects in the existence of a hidden symmetry for 
cyclic heat engines. That is, under the parity-time (V'T) transformation, the heat 
exchanges between the engine and its two reservoirs are dual. Here the parity inversion 
indicates interchanging parameters related to the hot reservoir (the quantities with 
subscript h) and those related to the cold reservoir (the quantities with subscript c). 
The time-reversal operation changes the sign of time. This hidden symmetry can be 
mathematically expressed as 

VTQh = Qc, ^^<^'PTQc = Qh- ( 12 ) 

The above P7~-symmetry will be confirmed by the typical model of cyclic heat engines 
shown in the Sec. 13.31 and Sec. 13.41 It is not hard to check this point according to the 
following equations (l26|) and fl34|) . 

3.2. Constitutive relation for nonlinear response 

We consider a cyclic heat engine undergoes a thermodynamic cycle consisting of two 
“isothermal” and two adiabatic processes. The word “isothermal” merely indicates that 
the heat engine is in contact with a heat reservoir at constant temperature. In the 
process of “isothermal” expansion during time interval th, the engine absorbs heat Qh 
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from the hot reservoir at temperature T^. The variation of entropy in this process is 
denoted as AS. On the contrary, in the process of “isothermal” compression during 
time interval tc, the engine releases heat Qc into the cold reservoir at temperature Tc. 
There is no heat exchange and entropy production in two adiabatic processes. Assume 
that the time for completing the adiabatic processes is negligible relative to tc and th- 
So the period of the whole cycle is to = A + U- The heat exchanges Qh and Qc can be 
expressed as 


Qh = Tn{AS-S]:), and - Qc = Tc{-AS - (13) 


where (or S'*”) represents the irreversible entropy production in the “isothermal” 
expansion (or compression) process. This model of cyclic heat engines is of broad 
generality if the entropy production in adiabatic processes of the thermodynamic cycle 
can be neglected. Equation flT^ shows that Qh and Qc indeed abide by the same 
function type. 

The “isothermal” expansion (or compression) process may be regarded as reversible 
in the long-time limit th ^ oo (or tc —>• cxd). That means, Sh (or S'*”) should be vanishing 
when —)■ oo (or tc —)■ oo). If the entropy production is an analytical function, we may 
write out 


oiT 


^+Ah 

th 




and Sf 





, (14) 


with the time-independent parameters T/i, Ah, Tc, and Ac. The parameter Ah either 
depends merely on the detailed “isothermal” expansion process, or equals to some 
constant independent of specific processes. Similarly, Ac either depends merely on the 
detailed “isothermal” compression process, or equals to some constant independent of 
specihc processes. Up to the hrst order, equation (ITT)) degenerates into the so-called 
low-dissipation assumption proposed in [33] . 

Substituting (1TT|) into flT^ . we obtain the expressions of the heat exchanges 


Qh = ThAS - ThTh/to - ThTlAh/tl + 0{l/tl), 

Qc = TcAS + Tcfc/to + TcflAc/tl + 0{l/tl), 


with parameters f/j = ThH/th and fc = Tcto/A- 

The PT-symmetry 0121) between Qh and Qc, will impose a constraint on the 
coefficients Ah and Ac in ([15]). Substituting flT^ into flT^ . we hnd the only choice 
is 


A;, = -Ac = -A, (16) 

where A is a process-independent parameter, and of course independent of the weighted 
parameters Sh and Sc which represent the degree of asymmetry of interactions between 
the engine and the reservoirs. With ffT6l) . the heat fluxes between the engine and two 
reservoirs may be expressed as 

Qh = Qh/H = ThAS/to - ThTh/tl + 0{l/tl), 

Qc = Qc/to = TcAS/to + TcTc/tl + 0{l/tl). 


( 17 ) 
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Now, we will construct the mapping from this cyclic heat engine into the generic 
model for tight-coupling heat engines mentioned in Sec. [2j Considering flT7|l . the 
weighted thermal flux ([2]) may be further expressed as Jt = {scTh -|- ShTc)AS/to — 
{scThTh — ShTcTc)/tl + 0{l/tl). According to the physical meaning of Jt, which has 
been fully discussed in [36], the quadratic order term of 1/fo should be vanishing. Thus, 
we require ScThTh — ShTcTc = 0. Combining with Sh + Sc = 1, we obtain the weighted 
parameters 


Sh 


ThTh 

ThTh + TcTc 


and Sc 


TcTc 

ThT' h + Tc^c 


(18) 


With these weighted parameters and dehnition ([6|) for cyclic heat engines, we further 
derive the weighted thermal flux ([2|) and the weighted reciprocal of temperature ([3|) as 




(19) 


and 


(9 


h/i -|- Tc 

Thfh + Tcfc^ 


( 20 ) 


with the coupling strength ^ = ThTcfdAS. In this work, we focus on the constitutive 
relation accurate up to the quadratic order. Thus the third and higher order terms of 
Jm could be neglected. Within this scope, equation flT^ implies that Jt is still tightly 
coupled with Jm- From Eqs. flTdl) . flTdl) . flTHll . (1201) and dehnition = 1/to , we obtain 
the expression of the generalized mechanical force as 


Xm = -{Th - Tc)(3AS + (f, + f,) - {TtSl - TcTl)(3KJl + 0(4). (21) 

From the above equation we hnally solve the generic nonlinear constitutive relation for 
tight-coupling cyclic heat engines 

Jm = p p ^ [1 + A{sh - Sc)A] 0{A^, X^), (22) 

t /x + t c 

with the consideration of flT0|) . flTSl) . fl20l) and ^ = ThTcfdAS. 

For the low-dissipation heat engine [33], the entropy production flTT|) in the 
“isothermal” expansion (or “isothermal” compression) process merely contains the hrst 
order term of Th/th (or Fc/tc)- Thus, the parameters Ah, Ac as well as A vanish 
in this model. Then, equation fl22l) degenerates into a linear constitutive relation, 
Jm = A/(Th + Fc), which is consistent with the result derived in [36] . 


3.3. Curzon-Ahlhorn endoreversible heat engine 

The Curzon-Ahlborn endoreversible heat engine [3] undergoes a thermodynamic cycle 
consisting of two isothermal processes and two adiabatic processes. In the isothermal 
expansion process, the working substance is in contact with a hot reservoir at 
temperature Th. Its effective temperature is assumed to be The (The < Th). During 
time interval th, an amount of heat Qh is transferred from the hot reservoir to the 
working substance according to the heat transfer law 


Qh — l^h(Th — The)th, 


(23) 
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where Kh is the thermal conductivity in this process. The variation of entropy in this 
process is denoted by AS. In the isothermal compression process, the working substance 
is in contact with a cold reservoir at temperature Tc. Its effective temperature is Tee 
{Tee > Te). During time interval te, an amount of heat Qe is transmitted from the 
working substance into the cold reservoir according to the heat transfer law 


-Qc = l^c{Te - Tee)te, (24) 

where Ke denotes the thermal conductivity in this process. Expressions fl23|) and fl2T|) 
show that the laws of heat absorbed and that released by the engine indeed abide by 
the same function type. The heat exchanges and the entropy productions are vanishing 
in the two adiabatic processes. Besides, the durations of these two adiabatic process are 
assumed to be negligible comparing to te and th- Thus the period (to) for completing 
the whole cycle is te + th- 

According to the discussion in ca, the endoreversible assumption proposed by 
Curzon and Ahlborn [3] leads to 


AA = Qh/The = Qc/Tee. 

Thus, equations fl2^ and fl2Tl) may be transformed into 

ThAS 


Qh — 


-, and Qe = 


TeAS 


(26) 


(26) 


(27) 


l + AS/7»i(,’ 1-AS/7A’ 

respectively, with two parameters 'jh = i^hth/to and 7 c = Kete/to. It is easy to verify 
that the heat exchanges [equation fl26l) ] indeed satisfy the "PT-symmetry flT^ for cyclic 
heat engines. 

Now, we will map this Curzon-Ahlborn heat engine into the generic model detailed 
in Sec. |2]and derive the constitutive relation for nonlinear response. According to fl26ll . 
the heat fluxes may be expressed into 

Qh = Qh/to = ThAS/to - ThASQ'yttl + 0(1/4), 

Qe = Qc/to = TeAS/to + TeASQ^etl + 0(1/4). 

From ([2]), fl271) and the dehnition = 1/to for cyclic heat engines, we can derive the 

weighted thermal flux Jt = (seTh + ShTe)ASJm - (seTh/lh - ShTe/'ye)AS‘^+ 0(Jl/). 

Considering the physical meaning of Jt discussed in [36], we require SeTh/'-fh — ShTe/'jc = 

0. Combining with Sh + Se = 1, we derive the weighted parameters 

Thlc 1 Te'^h /nQ\ 

Sh = -——, and Sc = -—;—. (28) 

Thic + Te'^h ThJe + Te'Jh 

Then the weighted reciprocal of temperature ([3]) and weighted thermal flux ([2]) may be 
expressed as 

Jh + 7c 


/? = 


T/j7c + Te'Jh ’ 


(29) 


and 


J, = + 0(J’), 


(30) 
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with ^ = T/jTc/SAS'. As mentioned above, we just focus on the constitutive relation 
accurate up to the quadratic order. In this sense, Jt is still tightly coupled with 
Substituting (l26ll and (12^ into ([6]), we obtain the generalized mechanical force 

Xr^= - (T, - T,)I3XS + (T,/7, + 

- + 0( J^), (31) 

from which we hnally solved the constitutive relation for the Curzon-Ahlborn heat engine 
accurate up to the quadratic order 

^ (7e + 7h)A^2^ ^ ^ 

with the consideration of flTOD . (1^ . fl2^ and ^ = ThT^dAS. It is obvious that the 
constitutive relation fl3^ for the Curzon-Alhborn heat engine is a specihc form of the 
generic nonlinear constitutive relation fl22l) for tight-coupling cyclic heat engines with 
model-dependent parameters A = l/AS", Th = AS'^/'jh, and fc = AS'^/jc- 


3.4- Revised Curzon-Ahlborn heat engine 


The thermodynamic processes and dehnitions of physical quantities in the revised 
Curzon-Ahlborn endoreversible heat engine are exactly the same as those in the 
original Curzon-Ahlborn heat engine mentioned in Sec. 13.31 except for the heat transfer 
law in two isothermal processes. Here the law of heat exchanges is revised to 

Qh = i^hiTCe - TC^)th, and -Qc = - TC^)tc- (33) 


The above equation shows that the heat exchanges in two isothermal processes still 
abide by the same function type. 

Considering the endoreversible assumption (|25|1 . we transform (|3^ into 


Qh 


Hill 1 

2H 7i«o j 


and Qc 


Ictp 

2T, 




4A^T2\ 

IcH ) ’ 


(34) 


respectively, with the parameters jh = i^hth/to and 7 c = KcA/^o- If is easy to verify that 
the heat exchanges [equation fl341) ] indeed satisfy the PT-symmetry flT^ for cyclic heat 
engines. 

This revised Curzon-Alhborn heat engine may be mapped into the generic model. 
The mapping procedure is similar to that in Sec. 13.31 Then we obtain 


^ The ^ Thh 

TH + Thh’ “ Th, + Th,: 


(35) 


The + Thh 
T'he + ThV 


(36) 


J, = T^TJASJ^ + 0(4) = ^ J™ + 0(7„), 


(37) 


and the nonlinear constitutive relation accurate up to the quadratic order 


Ad. 


Ihlc 


(T^yc + T^^jQAS'- 


:A 


1 + 


2(sh - gc) 

AS 


A 


+ 0(A^Xi). 


(38) 
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Obviously, this nonlinear constitutive relation for the revised Curzon-Ahlborn heat 
engine is a special form of the generic constitutive relation fl2^ for tight-coupling 
cyclic heat engines with model-dependent parameters A = 2/AS, Th = T^AS"^/'jh, 
and r, = 

4. Autonomous heat engines 

In this section, we will investigate the hidden symmetry in autonomous heat engines 
and its influence on the constitutive relation for nonlinear response. 

4-1. V-symmetry for autonomous heat engines 

By analyzing typical models of autonomous heat engines in the literature such as the 
Feynman ratchet as a heat engine [SI] and the single-level quantum dot heat engine |32] , 
we found that the expressions of the forward and backward flows in an autonomous heat 
engine conform to the same function type. For example, the forward and backward flows 
in the Feynman ratchet as a heat engine [31] may be expressed as Rp = and 

Rb = [36], respectively, where the detailed meanings of physical quantities 

will be fully explained in Sec. 14.31 Obviously, the forward and backward flows conform 
to the same function type. This homotypy reflects in the existence of a hidden symmetry 
under the parity {V) inversion discussed as below. The parity inversion indicates 
interchanging the quantities related to the hot reservoir (the quantities with subscript 
h) and the the quantities related to the cold reservoir (the quantities with subscript c), 
and simultaneously reversing the sign of the external load. Under the parity inversion, 
the expression of forward flow would turn into the expression of backward flow, and vice 
versa. This P-symmetry may be mathematically expressed as 

RRp = Rpi and RRpj = Rp. (39) 

On the other hand, the forward flow in an autonomous heat engine should merely 
rely on the intrinsic quantities related to the hot reservoir, such as Th, the temperature 
of the hot reservoir, and Sh, the relative strength of interaction between the engine and 
the hot reservoir. The backward flow should merely rely on the intrinsic quantities 
related to the cold reservoir, such as Tc, the temperature of the cold reservoir, and Sc, 
the relative strength of interaction between the engine and the cold reservoir. Thus, 
combining with the 7^-symmetry fl3^ . the forward flow Rp and backward flow Rp rnay 
be formally expressed as 

Rp = ^{shXm, I/Th), and Rp = ^{-ScXm, I/UJ, (40) 

respectively, where <l>(a:, y) represents a function type with independent variables x and 
y. The opposite signs before the terms related to Xm in Rp and Rp originate from the 
fact that the external load changes its sign under the parity inversion. The P-symmetry 
fl3^ as well as the formal expression fl40|) of the flows will be conhrmed by the examples 
of autonomous heat engines shown in Sec. 14.31 and Sec. 14.41 
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Figure 2: (a) Schematic of a generic autonomous heat engine; (b) Schematic of 

reduced autonomous heat engine. 


The schematic of a generic autonomous heat engine is depicted in Fig. [2]^a). In 
this setup, the heat engine consist of + 1 separate states Sn (n = 1, 2, ■ • •, iV + 1), in 
which tSi denotes the hot reservoir at temperature Ti = and iStv+i denotes the cold 
reservoir at temperature T/v+i = T^. The particles are transmitted sequentially from 
the hot reservoir tSi into the cold reservoir S^+i through the iV — 1 intermediate states. 
Thus the whole transmission process TZ contains N elementary transmission processes 
7^* (i = 1, 2, ■ • •, N) between the adjacent states as shown in Fig. [2]^a). Each elementary 
process TZ^ is composed by a forward process with flow Rp and a backward process 
with flow Rp. It should be noted that Rp (or Rp) is the product of the pure forward 
(or backward) jumping rate in elementary process TZ^ and the occupation probability of 
state Si (or tSj+i). When the engine operates in the steady state, the relation 

Wp - R^p = r, (41) 

is satished in each elementary process TZ^ {i = 1, 2, • • ■ , N), where r is the net flow of 
the whole transmission process TZ. In this engine, along the particle flow r, the heat 
absorbed from the hot reservoir iSi is partly released into the cold reservoir and 

partly outputted in the form of work against a globally generalized external force Xm- 
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It has been proved that the entropy production rate may be expressed as a = 
rJ2f=i^^{RF/R b) autonomous heat engines operating in a steady state 

Comparing it with (|TT|1 . we obtain the expression of affinity 


N TDi 

A = Vln^, 

i=l 


(42) 


by considering dehnition = r for autonomous heat engines. Substituting fHTD into 
(02]), the affinity may be further expressed in terms of or as 


R\ 


N 

A = E In 

i=l ' 

^ _|_7’ 

A=Eln^ 

1=1 ^ 


N N 

= E -^r + I E 


■ 1 Rp 

l=l 2- 


2 ^ Ri 

t=l Ftp 


N N 

= E -^r - I E 


• 1 Rr 

1=1 F 


2 ft' 

1=1 Ftp 


rr^ + O(r^), 
+ 0{r^), 


(43) 


respectively, where O(r^) denotes the third and higher order terms of r. 

Assume that the forward and backward flows in all elementary transmission 
processes 72.* {i = 1,2, ■ ■ ■, N) conform to the same function type. Focus on one 
elementary transmission process 72.* between state Si at effective temperature Tj and 
state iSj+i at effective temperature E+i, which could be regarded as a minimal setup 
of autonomous engine operating between “hot” reservoir Si and “cold” reservoir iSj+i. 
According to f0O|) . 1/72^ and 1/7?^ can be formally expressed as 

l/Rp = l/E) = (j)i + (l)'iS\Xl,^ + 0{XI^ ), 

^/Rb = l/^*+i) = </>*+! - ^'i+A^\n + 0{X\^), 

respectively, where s\ and s* (satisfying + s* = 1) are weighted parameters in this 
minimal setup of heat engine, which denote the asymmetry degree of interaction in 
elementary transmission process 72.*. X\^ satisfying Eili denotes the locally 

generalized external force in process 72.*. Then we can dehne the weighted parameters 
for the whole autonomous heat engine as 


Sh = 




N 


and Sr = 




N 

■ V s*X* 


‘-m i=i i=i 

In (041) we have dehned the parameters and 0' as 
1 

(x=0,y=l/Ti) 


$ 


and 0' = 


dx 


{x=Q,y=l/Ti) 

If we further introduce an average quantity of (/>' (z = 1, 2, • • ■, X + 1) as 


1 N+l 

N+ik 


it is easy to verify 


and 


b' = 0{Xt) (1 < z < X + 1), 


(45) 


(46) 


(47) 


(48) 


4>i - 4’i ^ 0{Xt) {l<iJ<N + l), 
from definitions and 1/T/. < l/T, < l/T^. 


( 49 ) 
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Substituting 


^ = b E In 


into (H5]) . we obtain the expression of A as 


\i=l 

N ± 


— r 


= E 


i-p 

+ 0i+l 


+ ^ln. 


R\ 


i=l 

U 


'vi=l 


'■+ - Sc)X„r + 0{X^,r). 


(60) 


From fl50p and definition = r, we can solve the constitutive relation for nonlinear 
response 


Jm. — LA 


■("S/i Sc)X^ 


+ 0{A\X-^), 


(51) 


with the parameter L = + 0i+i)/2] 


4-3. Feynman ratchet 


r 



Figure 3: Schematic digram of Feynman ratchet as a heat engine. 

The Feynman ratchet [TTI - HS] may be regarded as a Brownian particle walking in a 
periodic potential with a hxed step size 9. As depicted in Fig. El the Brownian particle 
is in contact with a hot reservoir at temperature Th in the right side of each energy 
barrier while it is in contact with a cold reservoir at temperature Tc in the left side of 
each barrier. The particle moves across each barrier from right to left and outputs work 
against an external load 2 ;. The height of energy barrier is e. The width of potential in 
the left or right side of the barrier is denoted by 9c or 9^ = 9 — 9c, respectively. In the 
steady state and under the overdamping condition, according to the Arrhenius law, the 
forward and backward flows can be respectively expressed as 

Rp = and Rb = (52) 

where ro represents the bare rate constant with dimension of time“^. According to fl52ll . 
the forward and backward flows in the Feynman ratchet indeed conform to the same 
function type. The net flow in the Feynman ratchet may be dehned as r = Rp — Rb- 
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In each forward step, the particle absorbs heat e + z6h from the hot reservoir and 
releases heat e — z6c into the cold reservoir, while outputs work w = z9 against the 
external load. The energy conversion in each backward step is exactly opposite to that in 
the forward step. Thus, the heat absorbed from the hot reservoir and that released into 
the cold reservoir per unit time, as well as the power output may be further expressed 
as 

Qh = {e-\-z9h)r = er + {eh/9)W, , . 

4 = (e - z9,)r = 6r - {9J9)W, ^ ’ 

and 


W = z9r, (54) 

respectively. Comparing fl53l) with (jl]), we can straightforwardly write out the weighted 
parameters 


Sh = 9hl9, and Sc = 9cl9. (55) 

Substituting fl5^ into ([2]) and ([3]), we derive 

Jt = er, and /? = {9h/Th + 9jT,)/9, (56) 

from which we derive .^ = e by using ([9]) and dehnition Jm = ^ for autonomous heat 
engines. From ([7]) and fl541) . we obtain the expression of generalized mechanical force 


= -I3z9. 

Considering (l5^ -(l57 )) , we transform (l52|) into 


(57) 

(58) 


It is easy to verify that the forward and backward flows in fl58l) indeed satisfy the 
P-symmetry fl39l) and conform to the formal expression fl40|) . 

Substituting (ED, (HUD , (EZD and (|58ll into Jm = r = Rp — Rp, we obtain the 
nonlinear constitutive relation up to the quadratic order as 


Jm = roe 


1 “1“ 2 ^c)^r7 


+ 0(A^Xi), 


(69) 


where (d = (l/T/j + Obviously, equation flSUP is a specihc form of the generic 

nonlinear constitutive relation flUTll for autonomous heat engines. 

From another point of view, we can also derive the constitutive relation fEUP from 
651. Comparing the forward and backward flows of the Feynman ratchet [in equation 
(EHD] with the formal expression (HUjl of forward and backward flows for autonomous 
heat engines, we could derive the specihc form of function ^{x,y) for the Feynman 
ratchet as 


y) = roe™'''’-'". 


(60) 


Because the Feynman ratchet may be regarded as a setup consisting of only one 
elementary transmission process between two states (i.e., two reservoirs), we use 
subscript h to denote the quantities related to the hot reservoir and subscript c to 
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denote the quantities related to the cold reservoir. Then, substituting fEU]) into 
and flT7|) . we obtain 


(t>h = — 

^^0 


= 

ro 


and 


0' = + €) = 


2rof3 \Th 




(61) 


(62) 


Finally, substituting fl6ip and fl62p into the generic nonlinear constitutive relation (IHip 
for autonomous heat engines, we can achieve the nonlinear constitutive relation fl5^ for 
the Feynman ratchet again. 


4-4- Single-level quantum dot heat engine 

A single-level quantum dot heat engine [32] consists of three parts: a hot lead at 
temperature and chemical potential Hh] a cold lead at temperature T^. (Tc < T^) 
and chemical potential fic (/^c > h'h)] a single-level quantum dot with energy level e 
(e > /ic), which located between the two leads. In the forward process, an electron 
jumps from the hot lead to the cold one via the quantum dot. The electron absorbs 
heat Qh = e — fih from the hot lead and releases heat qc = e — yic into the cold one, 
and simultaneously outputs work w = i^c — fih- In the backward process, the energy 
conversion exactly opposites to that in the forward process. When the engine operates 
in a steady state, the overall forward and backward electronic flows may be expressed 


a 


rp = 


and vb = 


a 


(63) 


respectively, with the coefficient a. The above equation indicates the overall forward 
and backward flows of this single-level quantum dot engine indeed conform to the same 
function type. Then the net flow from the hot reservoir into the cold one can be written 
as 

1 1 


Jm — B ^ 


(64) 


g(£ Uh)/Th _j_ ]_ g(e Uc)/Tc _j_ 

The heat absorbed from the hot reservoir and that released into the cold reservoir per 
unit time, as well as the power output could be further expressed as 

Qh (^ h'h^J’nm Qc (^ h‘c)Jrm ( 65 ) 

and 

W = (nc-iah)Jm, ( 66 ) 

respectively. As the engine operates in the steady state, the quantum dot is assumed 
to be locally in equilibrium. Thus we may introduce /i {qih < /x < /^c) as the effective 
chemical potential of the quantum dot. Then we can transform (l65ll into 

\ JTT /\ f \ J ( 


Qh 


Uc ^'h 


•} Qc 




ly. (67) 
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( 68 ) 

(69) 


Jt = iJm, (70) 

with the coupling strength ^ = e — ja. Substituting fl68|) into ([3]) and ([7]), we derive the 
expressions of weighted reciprocal of temperature and generalized mechanical force as 


/3 — Sh/Th + Sc/Tc — 


l^c f^h 


^ /^c 


n 


T. 


and 


a 


(71) 

(72) 
could be 

(,[ie-ti)-s,,Xm/l3]/n + 1 ’ ^ e[(s-ll)+S,X,r./l3]/Tc + 1 • 

Obviously the above equation satishes P-symmetry (15^ and abides by the formal 
expression (HOl) of the forward and backward flows. Besides, the specihc form of 
function ^{x,y) for single-level quantum dot engine can be written as ^{x,y) = 
aj Ue-r-)y-xy/y ^ x' _ 

Substituting flH5|l -fl72 |) and flTOl) into fl6^ . we hnally obtain the constitutive relation 
for nonlinear response as 

a 


Xm (dw (^f^h h')/7h (/^c h')/7c' 

With dehnitions fl68|l and fl72l) . the overall forward and backward flows in 
transformed into 

V TaTFF. i 7 1 Snd Vb 


Jm, 


4cosh (/3^/2) 


-A 


1+ 


+ 0(A^,Xi), (74) 


with fd = (T^ ^)/2. Obviously, the above equation is a specihc form of the nonlinear 

constitutive relation 051 p for tight-coupling autonomous heat engines. 


4-5. Reduction of autonomous heat engines 

In the steady state, the entropy production rate of the multi-step autonomous heat 
engine shown in Fig. Ei: a) may be expressed as cr = r hi(i?^/i?^), where r = 

Rp — R^p is satished in each elementary transmission process (i = 1, 2, • ■ •, N). From 
the perspective of entropy production, if we dehne an effective transmission process TZ 
with forward rate Rp and backward rate Rb, which satisfy 


R, 


Rb = r, and In In —^ = A, 

r? ^ R* 

fCb i=l FLb 


(75) 
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the multi-step model of autonomous heat engines shown in Fig. [2I^a) may be reduced to 
an effective model shown in Fig. Et^b) with effective forward flow Rp and backward flow 
Rb- Then, from dehnition (j75|l . we can derive the expressions of 1/Rf and 1/Rb as 


1 


Rf 


^ 1 
= y — 

r>i 


• 1 ^F 

2 = 1 ^ 


1 ^ 


1 


1 ^ 
-r + k E 


1 


N N 

J- = y -L 4-1 Y — 4-1 y 


® ^F^F^F 


Rb iYl ' 2 R'gR^g 




r^ + 0{r^), 
-|- O(r^). 


(76) 


From flT2|) and ([75]), the affinity may be expressed as A = ln[RF/ {Rf — r)] = 
r/Rp + r'^/2R\ -|- O(r^) or A = ln[(i?s -|- r)/RB] = r/Rp — r^/2R^ -|- O(r^). From 
these expressions, we may derive the affinity and entropy production rate as 


and 


A’’ 

2 Rp-r 


O' — Jrn.A 


1 Rp + r 1/1 1 \ 3 n 

- hi —^r + O(r^), 


R 


B 


2 \Rf 


R 


B, 


R] 4 + 0(4), 


(77) 


(78) 


2 Vflf 

respectively. Surprisingly, from (ITT]) and (176]) . we can prove l/Rp-hl/Rp = I^{^i(1/Rf+ 
1/7?^) -|- O(r^). Thus fITT]) and fITH]) may be further transformed into 


N 




2 S V«: 



and 



(79) 


(80) 


respectively. 

As an example, we will give the reducing procedure from the single-level quantum 
dot engine depicted in Sec. 14.41 which contains two elementary processes (7Z^ and 7Z^), 
into a reduced model containing only one effective process (7Z). As depicted in [32] . 
when the engine operates in the steady state, the forward and backward flows in the 
process between the hot reservoir and the quantum dot as well as those in the process 
between the quantum dot and the cold reservoir could be expressed as 


and 


ah _ ^h(l - fh) + ac(l - fc)^ f 

~ -^- (^hjh-, 

O^h + flc 


ah _ 

Kb — 


O^hfh T ®c/c r \_ 

' ®/i(l Jhjj 

O^h + 0-c 


(81) 




- f4 Ii% = ALJhUytLJAaJ^, (82) 


respectively, with fh = 1/[1 + and fc = 1/[1 + ah (or Oc) is a 

coefficient in process 77^ (or 77'^). 

Because there are only two elementary processes in the single-level quantum dot 
model, the second and higher order terms in 0761) is vanishing. Substituting 08ip and 












Hidden symmetries and nonlinear constitutive relations 


18 


fl52]) into (ESI), we derive the reciprocal of effective flows: 

1 ^ alfhjl - fh) + Qe/c(l - fc) + 2ahacfh{l - fc) - {ah + ac)r 

Rp O'hO^ci^Oihfh R (^cfc)\,^ {^(^hfh “1“ (^cfc)/(^(^h R ®c)]/'/i(l fc) 

1 ^ Qfe/fe(l - fh) + Qe/c(l - fc) + 2afeaJc(l - fh) + {ah + ac)r 

-Rg “1“ ®c/c) [l {^(^hfh “1“ Ojcfc) / iflh ®c)]/'c(l //i) 

Then, with consideration of a = ahad{au + Oc) and r = a{fh — fc), we conld verify 
Ri? — Rb = r and a = {Rp — Rp) iTiRp/Rp are exactly satisfied in this rednced model. 
Of course, both (1771) and fl78|) still hold. 

Further, with the consideration of ([S]), fIBSD . (ITT]) and (1771) . fh = 1/[1 + 
and fc = 1/[1 + may be expressed as 

^ 1 + f,i^-shXm/y){y-xt/2) ’ = 1 + (,{i+s.Xm/y){y+xt/2) ’ 


respectively, with jd = {1/Th + llTf)l‘l and ^ = e — [i. Substituting (IMll into (18^ . 
through tedious calculations we derive 


1 1 
Rp Rp 


-cosh^(/3^/2) 

a 



tanh(/5^/2) 

2 


(S/l Sc)Xyyi 


+ 0{XdX/), 


(85) 


with consideration oi fd/fd = 1 + 0(Xt) and ahac/{ah + Oc) = a. Then, combining 
(ITT)) and fl85p . we can achieve the nonlinear constitutive relation (1711) for the single-level 
quantum dot heat engine again. 


5. Conclusion 

In this work, we first investigated the hidden symmetries existing in a broad class of 
heat engines. In cyclic heat engines, this hidden symmetry may be characterized as the 
duality of the heat exchanges between the engine and two reservoirs under the parity¬ 
time transformation [RT-symmetry (IT^ ]. Based on a generic model of nonequilibrium 
cyclic heat engine and this RT-symmetry we derived generic nonlinear constitutive 
relation (12^ for cyclic heat engines. In autonomous heat engines, the hidden symmetry 
may be characterized as the duality of the forward and backward flows under parity 
inversion [R-symmetry (l39l) ]. which leads to formal expression (11(11) of the forward and 
backward flows. By applying this formal expression in a generic multi-step autonomous 
heat engine, we derived nonlinear constitutive relation (|5T1) for autonomous heat engines. 
The hidden symmetries as well as the nonlinear constitutive relations are all confirmed by 
typical heat engines in the literature that we have known. Besides, we also proved that 
the multi-step autonomous heat engine could be reduced to an effective autonomous 
engine containing only single transmission process from the perspective of entropy 
production. 
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